How wing compliance drives the efficiency of self-propelled flapping flyers 
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Wing flexibility governs the flying performance of flapping wing flyers. Here we use a self-propelled 
flapping-wing model mounted on a "merry-go-round" to investigate the effect of wing compliance 
on the propulsive efficiency of the system. Our measurements show that the elastic nature of the 
wings can lead not only to a substantial reduction of the consumed power, but also to an increment 
of the propulsive force. A scaling analysis using a flexible plate model for the wings points out that, 
for flapping flyers in air, the time-dependent shape of the elastic bending wing is governed by the 
wing inertia. Based on this prediction, we define the ratio of the inertial forces deforming the wing 
to the elastic restoring force that limits the deformation as the elasto-inertial number M&%. Our 
measurements with the self-propelled model confirm that it is the appropriate structural parameter 
to describe flapping flyers with flexible- wings. 



Flapping flight is probably the way of locomotion using 
the most complex dynamics in the animal realm [l], 
Structural properties of animal wings, together with wing 
kinematics, constitute the basic elements of a tough prob- 
lem: being able to stop, accelerate, execute sharp turn- 
ing, hover, etc... However, flapping costs a large amount 
of energy due to the perpetual cycle of acceleration and 
deceleration involved in the process of generating useful 
aerodynamic forces |3j. As a means of minimizing this 
cost, natural systems have presumably optimized animal 
wings by tuning their flexibility, which enhances not only 
their mechanical resistance, but also the animals' flight 
efficiency. The crucial nature of the elastic response of 
the wings in the propulsive performance of a flapping 
flyer has been made clear not only by observing nat- 
ural systems [3], but also investigating simplified mod- 
els where wing compliance determines drastic changes in 
thrust production and efficiency (see [4-8]). A few ways 
by which wing flexibility is favorable for both flying ani- 
mals and man-made devices have recently been proposed 
(9I-HH (see also extensive review by Shyy et al [12]). 
However, although a common hand waving argument is 
that wing compliance can be beneficial for the flapping 
flyer if elastic potential energy can be stored when the 
wings bend and released in a favorable part of the flap- 
ping cycle, the details of the balance of fluid dynamical, 
structural and inertial forces and moments that governs 
these mechanisms remain not well understood. Experi- 
mentally, one has to note that most studies do not con- 
sider self-propelled objects but the interaction between 
flapping bodies held static and an oncoming uniform flow 
that is driven independently of the flapping motion. If 
one thinks of cruising flapping flight, the very fact of de- 
coupling the flapping dynamics and the forward speed 
makes it difficult to extrapolate any conclusions about 
flight performance to the case of a free flying animal or 
machine. A notable exception is the experiment by Van- 
denberghe, Zhang and Childress [13], where a heaving 
wing mounted on an axis free to rotate was shown to 




FIG. 1: Experimental setup, (a) Sketch of the "merry-go- 
round". The wings are mounted perpendicularly to the ra- 
dial mast (length R — 0.5m). (b) Force measurement with 
a tabletop balance, (c) Close up view of the flapping de- 
vice, (d) Superposed pictures of the flexible flapping wing 
showing the chord- wise deformation. The semicircular poly- 
vinyl- chloride (PVC) wings had a c = 30mm maximum chord 
(at mid-span). A fiberglass structure attaches each wing to 
the flapping mechanism and rigidities the leading edge. The 
peak-to-peak flapping amplitude A at the span position of the 
maximum chord was set to A — 30mm. 

spontaneously give rise to a cruising speed perpendicular 
to the direction of the heaving motion. This work has 
been recently extended by introducing a pitching degree 
of freedom to mimic wing compliance |l4j . 

Here we use an experimental self-propelled flapping- 
wing model to study the effect of wing flexibility on its 
propulsive performance. The setup has been designed to 
enable measurements of the cruising speed, the thrust 
force, as well as the consumed power, as a function of 
the imposed wing motion (flapping frequency) and wing 
design. It is shown that increasing wing flexibility leads 
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FIG. 2: (a) Electrical power consumed by the motor (b) cruising velocity and (c) thrust force as a function of the napping 
frequency for the four tested pairs of wings. In (a) the points correspond to the force measurements with the system held at a 
fixed station, whereas the lines are the measurements taken with the system turning at its cruising speed. 



to a substantial reduction of the consumed power as well 
as an increase of the thrust power, both quantities being 
directly dependent on the flexural properties of the wings. 

The experimental setup consists of a flapping- wing de- 
vice that is allowed to turn on a "merry-go-round " type 
base (in the spirit of le petit manege from Marey [ 1 5l| ) . 
The two-wing flapper is attached to a mast that is ball- 
bearing mounted to a central shaft in such a way that 
the thrust force produced by the wings makes the flap- 
per turn around this shaft (see figure [J). For the chosen 
wing geometry (half disk of diameter S = 2L = 6 cm), 
the control parameters of the experiment are the flapping 
frequency (/) and the foil chord- wise flexibility (governed 
by its thickness h). Four pairs of wings were tested, of 
thicknesses 0.15, 0.23, 0.25 and 0.4 mm that correspond 
to masses per unit area of 0.20, 0.30, 0.33 and 0.53 
kg m -2 and bending rigidities (B) of 6, 25, 34 and 120 
mN m, in the range of the flexural stiffness of real insect 
wings (l6| . The natural frequencies of oscillation of the 
wings measured with relaxation tests are 71, 111, 125 
and 166 Hz, respectively. The measured quantities are 
the power consumption (P^, computed from voltage and 
current measurements in series [see also [l?]]) from which 
the power of the system running with no wings has been 
subtracted, the cruising speed of the device (U obtained 
using the measured time per revolution T = 2ttQ~ 1 ) 
and the thrust force Ft- The force measurement was 
performed by holding the flapper in a fixed station us- 
ing a string attached through a pulley to a calibrated 
weight as shown in figure [1] (top right) and monitoring 
the weight deficit on a tabletop balance as a function of 
the flapping frequency. The experiments reported here 
were performed with flapping frequencies ranging from 
10 to 30 Hz, which drove cruising speeds between 0.2 
and 1.5 m s _1 . Thus, the chord-based Reynolds num- 
ber Re = XJcjv (where c is the maximum chord and v is 
the kinematic viscosity) ranged between 1000 and 3000 
whereas the amplitude-based Strouhal number usually 
used to characterize flapping-based propulsive systems 



St a = fA/U ranged between 0.3 and 0.9 (using the flap- 
ping amplitude at mid-span to define A). These values of 
Re and St a are in the ranges that correspond to flapping 
flight in nature. 

The measurements are summarized in Fig. [2j First, 
the electrical power consumed by the motor Pi as a func- 
tion of the flapping frequency / is shown in Fig. [2] (a). 
As can be observed, keeping the system at a certain fre- 
quency requires less power as the wings become more 
flexible. It is also worth noting that the consumed power 
Pi as a function of the forcing frequency does not change 
whether the system is running (lines) or not (points), 
i.e. whether we are measuring the thrust force or the 
cruising speed. This suggests that the modification of 
the aerodynamic forces due the the cruising motion does 
not play a crucial role in the energy needed to perform 
the flapping. We note that the most flexible wing tested 
can save up to 60 % in consumed power with respect to 
the most rigid one in the 25 — 30Hz frequency range. 
Figures [2] (b) and (c) respectively show the cruising for- 
ward flight velocity U = Rfl and the thrust force Ft as a 
function of the input power. Again, for the wing rigidities 
tested here performance increases with increasing flexi- 
bility. In a general way, the present results show that 
the passive mechanisms associated with wing compliance 
can increase flight efficiency by increasing the thrust force 
and the cruising velocity while spending less energy. This 
efficiency enhancement is summarized in Fig. [3j where 
the ratio of thrust power P T = F T x U to input power Pi 
is displayed. All measurements are thus combined in this 
efficiency factor fj = Pt/P%, which has been normalized 
in the plots with respect to its maximum value for the 
most rigid wing as r] = fj/ max/frigid- Fig. [3] (a) shows 
that each wing has an optimum flapping frequency be- 
yond which the efficiency starts to decrease. It is worth 
noting that this optimum occurs significantly below the 
natural relaxation frequency so that any conclusion on 
the role of a resonance to minimize the cost of bend- 
ing (as discussed for instance in the case of undulatory 
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FIG. 3: Normalized efficiency factor 77 as a function of (a) fre- 
quency and (b) Strouhal number StA = fA/U. The efficiency 
factor 77 = Pt/P% has been normalized as r\ — fj/ max frigid- 



propulsion by [18]) is not applicable (see also [19|)- The 
same efficiency factor is plotted in Fig. [3](b) with respect 
to the Strouhal number StA, where now the position of 
the optimum decreases with wing flexibility. It occurs on 
a range between 0.3 and 0.45, which is consistent with 
the range of optimal Strouhal numbers observed in na- 
ture [20]. The 77 vs. StA curves show different behav- 
iors depending on the bending rigidity of the wing and 
prompt us to look for a new scaling including not only 
aerodynamical variables but also the wing structural pa- 
rameters. 

The first point to address when considering the effect 
of wing flexibility is to identify the main force that will 
bend the wing. In the dynamic regime, an elastic flap- 
ping wing is subjected to both the fluid dynamic pressure 
acting on the surface of the wing and the inertia force due 
to the oscillating acceleration. A measure of the impor- 
tance of these two bending forces can be given using a 
simplified model for the flapping wing as a plate of length 
L, mass surface density p s and bending rigidity B (for a 
plate of thickness h and Young's modulus E, B ~ Eh 3 ) 
whose leading edge is heaving sinusoidally with frequency 
uj and amplitude A. The moment of the mean fluid pres- 
sure force scales then as Mf ~ pfU 2 L 3 = pfUj 2 A 2 L 3 , 
where pf is the fluid density and Uf = Auo is the max- 
imum flapping velocity, whereas the moment of the in- 
ertia force scales as Mi ~ p s L 3 Auj 2 . The ratio of these 
two moments is actually a mass ratio -^j, which is 
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greater than 10 for all the wings tested in the present 
case. The main bending factor in this case is thus the in- 
ertia force, which will be counterbalanced by the elastic 
restoring force produced by the bent wing. This is con- 
sistent with the analysis by 21] who concluded for most 
wings moving in air that the feedback between fluid pres- 
sure stresses and the instantaneous shape of the wing is 
negligible with respect to the inertial-elastic mechanisms. 
The results in Fig. [2] (a) also support this assumption, 
since the change in the aerodynamic forces in the cruis- 
ing regime with respect to the fixed station operation of 
the system is undetectable in the consumed power vs fre- 
quency curves. We therefore proceeded to compare the 
moment of the inertial force Mi to that of the elastic 
restoring force that scales as M e ~ B. The ratio M^/M e , 
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FIG. 4: Schematic diagram of the redistribution of aerody- 
namic forces by a bending plate model (top) with respect to 
a rigid plate (bottom). Fr r is the thrust force and For the 
drag forces for a rigid wing, respectively. 
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FIG. 5: (a) Non-dimensional thrust power pr (and input 
power pi in the inset) and (b) efficiency factor as a function of 
elasto-inertial parameter Mei- The arrows in (b) indicate for 
each wing the value of M e % where the trend of the efficiency 
curve changes. The error in the Mei measurements is of 5% 
due to the uncertainty in the measurements of B. 



which we define as the elasto-inertial number 
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in terms of the bending length L5 = (B / PsAuj 2 ) 1 ^ 3 , mea- 
sures thus to what extent the inertial force due to the 
oscillating acceleration will be balanced by the elastic 
resistance to bending (analog definitions of this bend- 
ing length arise in problems where other forces drive the 
bending, see for instance |22l 123] for capillary and hy- 
drodynamical forces, respectively). This definition de- 
termines for instance that for M e % <C 1 the wing is too 
rigid for the inertia of the oscillating wing to have an ob- 
servable effect, or in terms of the bending length, L5 ^> L 
so that no deformation can be observed over the length 
scale L of the wing chord. 

Physically, the form of the bending wing can be seen 
as a "shape factor" that redistributes the contribution of 
the aerodynamic forces in both directions -of the flapping 
motion normal to the wings Fd [26| and of the forward 
displacement- as sketched in Fig HI During the flapping 
motion, the wings experience strong drag as they push 
fluid up and down during a stroke cycle. Because of the 
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flexibility of the wing, the experienced drag scales on a 
length depending on L& [231425], instead of L as it would 
for the rigid case. On the other hand, the change in shape 
induces a contribution of the aerodynamic pressure load 
in the forward direction that is also dependent on the 
wing bending. The two forces Fd and Ft, respectively 
normal and in the direction of the cruising speed, should 
therefore be directly dependent on the wing shape that is 
determined by the non-dimensional number M e i defined 
in Eq. [U This is clearly shown in Fig. [5] (a), where 
the useful thrust power and the input power are plotted 
as a function of M e i • The thrust power is defined by the 
product of the cruising speed and the thrust force (shown 
in Fig. [2] (b) and (c), respectively), which are rendered 
non-dimensional using the scalings, fr = FtL/B and 
u = U/Auj. The non-dimensional thrust power is thus 
defined as pt = UFtL/ BAuj. Fig. [5] (a) clearly shows 
that all the measured data collapse onto a single power 
law with J\f e i- The input power (shown in the inset of 
Fig. [5] (a)) is non-dimensionalized using the same scaling 
Pi = PiL/BAuj, and gives a measure of the work of the 
drag force experienced by the flapping (i.e. ~ Fjju^). 
This quantity also scales with the elasto-inertial number 
M e i- The efficiency (plotted in Fig. [5] (b) as a function 
of Mei) that we have defined previously as Pt / Pi is thus 
a ratio of works: the work of the useful thrust force Ft 
over the work of the drag force experienced by the flap- 
ping wing Fd- A salient feature of the 77 vs. M e % plots 
in figure [5] (b) is that the main trend followed by the 
experimental points over different sets of wings is valid 
only up to a certain threshold for each wing after which 
the measured data show a change of regime. We have 
noted in figure [3] that a simple resonance at the relax- 
ation frequency cannot explain the observed behavior. 
The possible role of subharmonic resonances in the ef- 
ficiency curves involving a detailed study of the phase 
dynamics is the subject of ongoing work. 

In summary, the effect of wing flexibility on the ef- 
ficiency of flapping flyers can be thought of as a two- 
step process: a solid mechanics problem, where the bal- 
ance between inertial and elastic forces determines the 
instantaneous shape of the flexible wings, followed by a 
fluid dynamics problem, where the boundary conditions 
set by the previous step govern the distribution of aero- 
dynamic forces. This simple passive mechanism, shown 
here to be well described using the elasto-inertial num- 
ber Mei, can bring a two-fold advantage: decreasing the 
energy cost while enhancing the thrust power. The self- 
propelled system described here gives a framework that 
should be useful to pursue further studies on the effect 
of structural and geometrical properties of the wings in 
the performance of flapping-based propulsion. 
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